Precision muon decay measurements and improved constraints on the weak interaction 
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The TWIST Collaboration has completed its measurement of the three muon decay parameters 
p, 8, and P M £. This paper describes our determination of p, which governs the shape of the overall 
momentum spectrum, and 8, which controls the momentum dependence of the parity-violating 
decay asymmetry. The results are p = 0.749 77 ± 0.000 12(stat.) ± 0.000 23(syst.) and 8 = 0.750 49 ± 
0.000 21(stat.)±0.000 27(syst.). These are consistent with the value of 3/4 given for both parameters 
in the standard model, and each is over a factor of 10 more precise than the measurements published 
prior to TWIST. Our final results on p, 8, and P M £ have been incorporated into a new global analysis 
of all available muon decay data, resulting in improved model-independent constraints on the possible 
weak interactions of right-handed particles. 

PACS numbers: 13.35.Bv, 14.60.Ef, 12.60.Cn 
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I. INTRODUCTION 



The TWIST experiment is a high-precision search for 
evidence of contributions to the charged-current weak in- 
teraction beyond those described by the standard model 
(SM) of particle physics. We take advantage of the purely 
lcptonic nature of the decay of the positive muon into a 
positron and two neutrinos, /i + — > e + v e 9^, which can 
be described to a good approximation as a four-fermion 
point interaction and in the SM is mediated by the W 
boson. 

The most general Lorentz- invariant, local, and lepton- 
number-conserving description is given by the matrix el- 
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7=S,V,T 
e,[i=L,R 
(n,m) 

where each scalar (S), vector (V), or tensor (T) inter- 
action between /x-handed muons and e-handed positrons 
has an associated coupling constant satisfying cer- 
tain normalizations and constraints [1]. Only 19 real and 
independent coupling constants are needed to describe 
entirely the interaction because g]^ R = and g\ L = 0, 
and a common phase is not observable. In the context of 
the V — A interaction of the SM, all coupling constants 
are zero except for g\ L = 1. The coupling constants pro- 
vide the probability Q eAt for a /x-handed muon to decay 
into an e-handed positron using 
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' e = fx and (J e „ = for e^//. In particu- 
lar, a model-independent limit on any muon right-handed 
couplings [l|, is determined from the probability 
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The differential muon decay spectrum [3(, using the no- 
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tation of Fetscher and Gerber Q , can be written as 
d 2 T 



P M cos 9 B F AS (x)} 



(4) 



where is the Fermi coupling constant, 9 S is the an- 
gle between the muon spin and the positron momentum, 
W efi « 52.8 MeV is the kinematic maximum positron 
energy, x — E e /W ef i is the positron's reduced energy, 
xq = rne/We/j. is the minimum possible value of x, cor- 
responding to a positron of mass m e at rest, and P M is 
the degree of muon polarization at the time of decay. P M 
is typically reduced from P", which is the helicity of the 
muon at the time of its production from a pion decay, 
due to depolarization undergone by the muon before it 
decays. 

The isotropic and anisotropic parts of the spectrum 
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are parametrized by four muon decay parameters p, r/, 
<5, and £, which are bilinear combinations of the cou- 
pling constants These four parameters, with the ad- 
dition of the radiative corrections Pjg C (x) and F^(x), 
are sufficient to describe the shape of the momentum- 
angle spectrum of the decay positron. We analyze the 
momentum- angle spectrum rather than the energy-angle 
spectrum out of convenience and because for these ener- 
gies the difference is insignificant. 

The introduction of chiral spin 1 fields to the SM has 
been investigated One consequence is that nonlocal 

tensor interactions appear, so that gJ^ L and g] iR are no 
longer zero. These new couplings can be measured in 
particular through the 6 parameter. 

Initial and intermediate measurements of p and 6 have 
already been published d-Q. This paper presents a de- 
tailed description of the final measurement of the p and S 
decay parameters by the TWIST Collaboration reported 
in Q . An identical and simultaneous analysis of the same 
data yielded the final P^t; parameter determination; a 
complete description with an emphasis on the system- 
atic uncertainties specific to P M £ was presented in [To| . 
The decay parameter r/ was fixed to the global analysis 
value of r\ — 0.0036 ± 0.0069 because the sensitivity 
to this parameter is reduced due to the multiplying factor 
x ~ 10~ 2 . 



II. EXPERIMENTAL SETUP 
A. TWIST spectrometer 

A brief description of the experimental setup is given 
here. A more detailed description of the apparatus can 
be found in [l2[ with the improvements made for this 
final analysis described in [Toj . 

An overview of the TWIST apparatus is shown in Fig. 
□ it was installed on the M13 beamline at TRIUMF, 
Vancouver, Canada. The 500 MeV proton beam from 
the TRIUMF cyclotron hit a carbon target producing pi- 
ons, some of which stopped and decayed near the surface 
of the target to create 29.79 MeV/c muons with 100% 
polarization. The beamline was tuned to transport these 
highly polarized muons with a central momentum of 29.6 
MeV/c and a momentum bite of 0.7% FWHM. The beam 
also contained several times as many positrons as muons, 
with the ratio varying with different tuning conditions. 
After passing through the beamline, the muons stopped 
at a rate between 2000 s" 1 and 5000 s" 1 in a thin target 
foil located in the center of the highly symmetric array 
of 44 planar drift chambers (DCs) [13| and 12 planar 
proportional chambers (PCs) composing the TWIST de- 
tector (Fig. [5]). The DCs and the PCs had an active 
region of 32 cm diameter and contained respectively 80 
and 160 parallel sense wires separated by 0.4 cm and 0.2 
cm. 
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FIG. 1. Conceptual drawing of the TWIST spectrometer. 
It shows the superconducting solenoid within the steel yoke, 
with the drift chambers and proportional chambers symmetri- 
cally placed upstream and downstream of the central stopping 
target. 

The DCs were filled with dimethyl ether gas and were 
assembled in modules of two or eight chambers in which 
the aluminized Mylar cathode foils were shared by neigh- 
boring chambers. DC 9-22 and 23-36, installed in two- 
chamber modules, formed a sparse stack covering most 
of the tracking region. The two eight-chamber DC 1-8 
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FIG. 2. (color online) A cross section of the TWIST detec- 
tor, including an example event of a downstream decay. PCs 
provide timing information and DCs determine the position 
of particles. The angle between the decay positron and the 
z axis, defined along the beam direction, is 6 — n — 8 S . The 
muon polarization direction is opposite to that of the z axis. 



and 37-44 modules instrumented the end of the tracking 
region [10]. PC 1-4 and 9-12 were installed at the ends 
of the detector for particle identification purposes. The 
PC 5-8 module had the target foil as central cathode foil 
and was installed in the center of the detector stack to 
make the entire array symmetric. The PCs were filled 
with a mixture of CF4 and isobutane. The array of low 
mass chambers was installed in a frame referred to as the 
cradle, filled with helium to further reduce the amount 
of material traversed by the muons and positrons. Two 
different target foils were used over two run periods to 
study the effects of the target material on the decay pa- 
rameters measurement: a (30.9 ± 0.6) /mi thick silver foil 
and the (71.6 ± 0.5) /jm thick aluminum foil used for the 
intermediate TWIST measurement [8j. Both metal tar- 
gets had purity exceeding 99.999% and featured minimal 
depolarization of the muons after stopping [l4j . 

The detector was installed in a superconducting 
solenoid producing a magnetic field of 2 T that was highly 
uniform over the tracking region and aligned with the 
beam direction. In order to obtain the required field uni- 
formity and also to reduce fringe fields, it was necessary 
to surround the solenoid cryostat with a cube-shaped 
yoke of approximately 3 m on a side. Two NMR probes 
were installed slightly beyond the radius of the tracking 
region in the cradle to monitor constantly the magnetic 
field strength during data taking. 

The magnetic field was mapped using a rotating arm 
equipped with Hall probes to measure the longitudinal 
component with a precision of 0.1 mT and an NMR 
probe for the total field. The Hall probes were sepa- 
rated by about 4.13 cm on the arm. A full rotation was 
performed every 5.0 cm along z for the central part of 
the tracking region, and every 2.5 cm for the edges of the 



region. The tracking region was fully mapped for each of 
the three field strengths used during data taking, 1.96 T, 
2.00 T, and 2.04 T. A smooth and higher granularity field 
map, including the relatively small transverse field com- 
ponents, was calculated using the Opera-3D software 
[l5j, matching the measured magnetic field map within 
±0.2 mT over the drift chamber region. 

The beamline vacuum pipe was extended through the 
fringe field region as close as possible to the end of the de- 
tector array. Upon exit from the vacuum, muons passed 
through elements of a "beam package," including a 20 cm 
length of gas degrader filled with an adjustable mixture 
of He and CO2 gas, a film strip degrader, and a muon 
scintillator that triggered the data acquisition system. 
The film strip degrader consisted of a roll of plastic film 
containing holes covered with Mylar degraders of varying 
thicknesses up to 0.1 cm. It could be rolled from outside 
the magnet yoke to choose which degrader was in the 
muon path. It was used to significantly degrade the muon 
beam momentum in order to stop muons well upstream 
of the target at the detector center, for special runs used 
for positron interaction studies (Sec. IIV A [I . The film 
degrader was set to an empty hole of the film strip for 
the normal acquisition of muon decay data. The muons 
traversed a total of ssl40 mg/cm 2 of material, including 
the beam package and the upstream half of the detec- 
tor, before stopping in the target. The transverse size of 
the beam spot was 1.6 cm FWHM. Because the cham- 
bers were operated at atmospheric pressure and thus the 
gas density varied with time, the ratio between the two 
gases in the gas degrader was automatically changed by 
a feedback loop to set and maintain the muon stopping 
distribution in the target. 

The downstream end of the detector was equipped with 
a second beam package during one data set to test the 
impact on the data of the asymmetry due to the presence 
of the upstream beam package. Two removable time ex- 
pansion chambers (TECs) were installed in the beam in 
the upstream fringe field region at the beginning and the 
end of each data set to characterize the muon beam prop- 
erties [16j . 



B. Experimental data 

The data used for the final phase were taken during 
fall 2006 for the Ag target and in summer 2007 for the 
Al target (see Tables Q] and HTJ the numbering of sets 
is not necessarily sequential). Monitor information was 
recorded during all runs for variables such as spectrom- 
eter temperatures, gas pressures and flows, and muon 
beamline element settings, and was later evaluated to 
identify any instabilities that could signify a low quality 
of data. Approximately 10% to 30% of runs in each set 
were discarded prior to the analysis to guarantee stable 
run conditions during the period of typically one week 
necessary to take a set. The criteria for rejection were 
conservative and unbiased; for example, they identified 
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runs with a problem in the data acquisition system, runs 
with a noisy chamber, or runs before the gas degradcr 
feedback loop was fully locked. 

The four nominal sets 74, 75, 84 and 87 were taken 
with optimal conditions for the measurement of decay 
parameters. For set 68, the degrader was changed so 
that the center of the muon stopping distribution was 
moved from near the middle of the target to a point only 
1/3 of the way through, to determine the sensitivity to 
stopping position variations. Set 83 was taken with a 
downstream beam package mirroring the upstream beam 
package to test the impact of the positrons backscattering 
into the spectrometer and the consistency of the results 
with or without a symmetric apparatus. Two sets (70 
and 71) were taken with different solenoid magnetic field 
strengths to verify that the decay parameters are insen- 
sitive to the transverse scale of the helices. 

Set 72 was unique in that it was taken with the TECs 
in place in the beamline, in order to test the effects of 
extra multiple scattering of the muon beam on the pa- 
rameter Pp£ through the depolarization of the muons, 
and also to monitor the stability of the muon beam posi- 
tion and angle over an entire week. The muon beam was 
steered off the detector axis with an angle y ~30 mrad 
for set 76 and with a position x w-1 cm and an angle 
6 X as- 10 mrad for set 86 to study the depolarization in 
the fringe field in simulation. Sets 70, 71, 72, 76 and 86 
were discarded from the P M £ measurement and used for 
systematic uncertainties studies due to their large depo- 
larization uncertainties [lOj ], but were used for p and 8 
since these parameters are insensitive to the muon polar- 
ization. 

The M13 central momentum was reduced to 28.75 
MeV/c for set 91 and to 28.85 MeV/c for sets 92-93 to 
study the effect of multiple scattering of the muons ex- 
iting the production target. The muons were stopped at 
the entrance of the detector for sets 73 and 80 by chang- 
ing the momentum selection and introducing a film de- 
grader in the beamline. These special sets of data are 
used to validate the simulation (Sec. IIV[) . 



III. ANALYSIS 

The muon decay parameters are extracted from the 
momentum-angle (p-0) spectrum of the decay positrons 
measured in the TWIST spectrometer. More precisely 
the difference in shape between the p-0 spectra from the 
data and from a full simulation of the TWIST apparatus 
is interpreted in terms of a difference in decay param- 
eters. A blind analysis is performed by using hidden 
decay parameters for the generation of the simulation 
17]- These parameters remain hidden until the end of 
the analysis when all systematic uncertainties and cor- 
rections have been determined to minimize the possibility 
that the results are affected by human bias. 

The simulation is analyzed using the same reconstruc- 
tion and event selection that is applied to the data, and 



TABLE I. List of Ag data sets used for the final TWIST mea- 
surement. The set numbers below are retained for historical 
reasons; missing numbers are not relevant for the analysis. 



Data Description 


Events (xlO 6 ) 


set 








oeiore 


Final 




cuts 


spectrum 


68 Bragg peak ^ 


1 1 1 
/ 41 


O A 


into target 






70 Central field at 1.96 T 


952 


50 


71 Central field at 2.04 T 


879 


45 


72 TECs in place, nominal beam 


926 


49 


73 Muons stopped at detector entrance 


1113 




74 Nominal 


580 


32 


75 Nominal 


834 


49 


76 Off-axis beam 


685 


39 


TABLE II. List of Al data sets used for the final TWIST mea- 


surement. The sets are listed in chronological order 


except for 


set 88, which was divided into the sets 91, 92 and 93 during 


the analysis because the running conditions changed. The 


set numbers below are retained for historical reasons; missing 


numbers are not relevant for the analysis. 






Data Description 


Events (xlO 6 ) 


set 








Before 


Final 




cuts 


spectrum 


80 Muons stopped at 


363 




detector entrance 






83 Downstream beam package in place 


943 


49 


84 Nominal 


1029 


43 


86 Off-axis beam 


1099 


58 


87 Nominal 


854 


45 


91 Lower momentum I 


225 


11 


(p = 28.75 MeV/c) 






92 Lower momentum II 


322 


15 


(p = 28.85 MeV/c) 






93 Lower momentum III 


503 


26 


(p = 28.85 MeV/c) 







reproduces very closely the detector response. Differ- 
ences between data and simulation arise from differences 
in the muon decay parameters and radiative corrections, 
and additionally from uncertainties in the simulation in- 
puts. The latter are the source of most of the systematic 
uncertainties. 



A. Simulation 

The Monte Carlo simulation of the TWIST experiment 
uses the geant 3.21 package jTs] to simulate the particle 
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interactions, the detector geometry, and its electronics. 
None of the physics processes undergone by the parti- 
cles such as bremsstrahlung or 5-electron production are 
modified or tuned from their definitions in GEANT 3.21. 
Since our apparatus had very thin scattering layers, for 
the energy loss we used the optional simulation of reduced 
Landau fluctuations with delta rays. 

The simulation includes all the elements necessary to 
reproduce accurately the muon and positron trajectories. 
The particles are transported in the Opera-3D magnetic 
field map using a classical fourth order Runge-Kutta nu- 
merical method. The description of the wire chambers 
includes the cathode planes and the wires, as well as their 
positions measured by the alignment calibration (Sec. 
IIIIF|) . The discontinuous behavior of the ionization of 
the wire chamber gas is simulated with ionization clus- 
ters generated randomly along the path of the charged 
particles. The ion cluster separation is matched to the 
data by comparing the timing of hits close to the wire 
in data and simulation. The drift time of each cluster is 
calculated from DC space-time relations (STRs) created 
by a Garfield simulation of the DCs. The effect 
of regions of the sense wires becoming temporarily in- 
efficient due to the presence of ionization from previous 
muon hits is also simulated. The data acquisition digiti- 
zation is part of the simulation in order to have output 
identical in format to that of the apparatus. 

For each data set, a corresponding simulation is gen- 
erated with its input parameters matched to the specific 
data taking conditions for that set, as needed. The frac- 
tions of He and CO2 in the gas degrader are set to time 
averaged values from the data. The muon beam pro- 
file measured by the TECs is used to generate the ini- 
tial muon directions [Ioj |. The muon and positron beam 
rates are matched to the data to simulate accurately the 
overlap in time of the hits in the DCs. Pions and cloud 
muonaj are beam particles that are not simulated be- 
cause they can be effectively eliminated from the exper- 
imental data. The magnetic field strength is matched to 
the cradle NMR probe measurements performed during 
each data set. Energy loss in some components outside 
of the tracking region is also simulated. For example, the 
upstream beam package had to be simulated in detail to 
reproduce the positrons scattering back into the detector 
and affecting the track reconstruction. The entire down- 
stream beam package was also included in the simulation 
matching set 83. 

Individual muons are generated at the location of the 
TECs, where the real beam has been well characterized, 
with polarization of 100% in a direction opposite to their 
momentum. The initial momentum and angle of the de- 
cay positron are generated with an independent program 



1 Cloud muons originate from pions decaying in flight as they move 
from the production target to the M13 beamline. These muons 
have a low polarization and are therefore removed during the 
analysis of the data with a time of flight cut. 



in order to isolate the hidden parameters of the blind 
analysis. The hidden parameters are chosen randomly 
within a range of ±10~ 2 from the SM values and remain 
encrypted during the whole analysis. The algorithm uses 
an accept-reject Monte Carlo technique with the theo- 
retical p-9 spectrum including full 0(a) radiative correc- 
tions with exact electron mass dependence, the leading 
logarithmic terms of 0(a 2 ), the next-to- leading logarith- 
mic terms of 0(a 2 ), leading logarithmic terms of 0(a 3 ), 
correction for s oft p airs and virtual pairs, and an ad-hoc 
exponentiation [2Cj]. The W boson's mass and the strong 
interaction contributions to the decay through loops are 
respectively on the order of 10~ 6 and 10 -7 [2l[, orders 
of magnitude smaller than our precision goal, and are 
therefore ignored for this measurement. 



B. Event and track reconstruction 

The reconstruction software is composed of three main 
algorithms. It begins by grouping the hits in the spec- 
trometer into different time windows and by identifying 
the type of particle (e.g., decay positron, beam positron, 
incident muon, secondary electron, etc.) causing the hits. 
Then a pattern recognition algorithm uses the positions 
of the hit wires to define helical tracks within each time 
window, using spatial information to separate the hits 
from two particles completely overlapped in time if nec- 
essary. Electron and positron tracks are finally recon- 
structed with high precision using the drift information 
in the DCs to extract the momentum and direction of 
the particles. 

Information from a 16 /xs interval is recorded for each 
event (from 6 /is before to 10 [is after the muon trig- 
ger) and divided into time windows designed to group to- 
gether the signals coming from each particle. The signals 
from the PCs define the beginning of the time windows 
because their time resolution is < 20 ns. The time win- 
dows are by default 1050 ns long to include the longest 
drift times in the DCs (50 ns before and 1000 ns after 
the first PC hit time). However if two particles are sep- 
arated in time by less than 1000 ns but more than 100 
ns, the first time window stops at the beginning of the 
second window. This type of event is rejected later in 
the analysis because signals of the particle in the first 
window can end up in the second window, confusing the 
track reconstruction. On the other hand, a time sepa- 
ration of less than 100 ns is not considered long enough 
for the PCs to identify two different particles and only 
one time window is created. In this case the signals cor- 
responding to each particle are separated by the pattern 
recognition using spatial information. This topology also 
includes the backscattcr of a decay positron from mate- 
rial outside the tracking region creating two independent 
tracks overlapping in time, as well as delta rays emitted 
in the tracking region. 

The particle identification algorithm uses the pulse 
widths in the PCs, roughly proportional to the energy 



6 



deposited, to separate muons from positrons since the 
two particles deposit different amounts of energy. Beam 
positrons are identified using the fact that they traverse 
the entire detector while the decay positrons originate 
from the target foil region in the middle of the chamber 
stack. The events are classified according to the particle 
content and the length of the time windows. 

The track reconstruction algorithm is performed on the 
signals in each time window. The first part of this algo- 
rithm is a pattern recognition, which combines hits on 
adjacent wires and associates signals together to form a 
coarse estimate of the helical track. The drift times are 
ignored at this stage and for this reason the Chebyshev 
norm is used as a fit optimizer [221 ] . This pattern recog- 
nition identifies and separates the tracks from the differ- 
ent particles contained in a time window, including J-ray 
electrons. A particle undergoing a large enough scatter- 
ing or energy loss due to the emission of a bremsstrahlung 
photon or a S-vay electron is reconstructed as two indi- 
vidual tracks by the algorithm. 

The next stage of the track reconstruction uses a x 2 
minimization to refine the helical trajectory identified by 
the pattern recognition. This helix fitter minimizes the 
residuals at each DC plane as well as kink angles in the 
center of each DC module, and includes as a fit param- 
eter the decay time of the muon. The time of flight of 
the decay positron to each DC plane is included in this 
calculation. The kink approach is well adapted to the 
TWIST spectrometer since the scattering masses are dis- 
crete [23| . The kink angles are weighted in the \ 2 min- 
imization by the inverse of the width of the Gaussian 
approximation calculated using the formula for multiple 
scattering through small angles [24| . For this analysis the 
space-time relationships used to convert the drift times 
into drift distances were measured using decay positron 
tracks (see Sec. MI Fl) . The trajectories between the DCs 
are calculated using the OPERA-3D magnetic field map 
to account for the inhomogeneities of the solenoid mag- 
netic field. The algorithm uses an arc step approxima- 
tion with variable size steps to integrate the magnetic 
field features. The energy lost by the positron through 
ionization is taken into account in the fitting procedure 
using 

cos ' 

i 

with AE the average energy loss of a track segment, k 
the thickness of the material i, and e*°™ the ionization en- 
ergy lost per unit of thickness in the material i calculated 
from the mean energy loss formulas [24j | . The track recon- 
struction has an inefficiency of a few 10 -4 , and an angle- 
dependent resolution at the end point (52.8 MeV/c), 
which is 58 keV/c when extrapolated to sin 6* = 1. From 
simulation, the absolute accuracy of the reconstructed 
momentum is better than 1 x 10 -4 . 



C. Event selection 

It is desirable to select classes of events that are very 
simple and therefore well simulated to reduce discrepan- 
cies between data and simulation. Our main selection is 
to find one muon and one decay positron separated by 
more than 1 (is. Events also containing a beam positron 
are kept only if the beam particle is separated from the 
incident muon and decay positron by more than 1 /is or 
less than 100 ns. A track from a decay positron backscat- 
tering at the upstream end of the detector and a beam 
positron track are indiscernible by the particle identifi- 
cation. The backscattering depends strongly on the de- 
cay positron momentum and angle. Thus events with 
a backscattered positron and events with overlap of de- 
cay and beam positrons within 100 ns are included in the 
analysis. These choices reduce the sensitivity of the anal- 
ysis to the accuracy in the simulation of these processes. 

The highly polarized surface muons are selected us- 
ing time of flight of the particles in the M13 beamline 
10]. A highly polarized muon beam is crucial for the 
measurement of P^, but also increases the sensitivity 
to the S parameter. The muons stopping in the target 
foil are selected by the next series of cuts. The first PC 
downstream and adjacent to the target acts as a veto for 
muons stopping too far downstream. The pulse widths 
in the two PCs just upstream of the target are used to 
eliminate muons that stopped in the gas or the wires of 
those chambers (IoT |. Also the muon position on the tar- 
get measured by the two PCs upstream is used to reject 
muons stopping more than 2.5 cm away from the central 
axis of the detector. Decay positrons from these rejected 
muons might not be contained within the tracking region. 

The purpose of the following selections is to identify 
which track corresponds to the decay positron. Tracks 
that failed the second stage of the track reconstruction 
and tracks corresponding to negatively charged parti- 
cles are rejected. The event classification determined on 
which side of the target the decay positron was emitted 
based on the side containing most of the hits. Tracks 
located on the opposite side are discarded. The next se- 
lection tries to match together tracks to check whether 
they originated from the same particle. In particular, the 
algorithm tries to match tracks from opposite sides of 
the target (using previously discarded tracks) to identify 
beam positron tracks and remove them from the analysis. 
In this case the criteria for a match are a time separation 
of less than 60 ns for the track times and a closest dis- 
tance of approach of the two extrapolations of the tracks 
of less than 0.5 cm. The matching can also identify tra- 
jectories split in two tracks (both located on one side of 
the target) due to a large scattering in a DC. In this case 
the closest distance of approach is only required to be 2 
cm. The position at the target of the muon as measured 
by the target PCs is compared to the extrapolation of 
the positron track back to the target to determine the 
vertex distance. An angle-dependent cut is applied to 
this vertex distance. If more than one track candidate 
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was selected, two more selections determine a single track 
corresponding to the decay positron. The tracks that are 
farthest from the target plane are discarded. If multi- 
ple tracks are equally close to target, the selected track 
candidate is the one with the shortest muon-positron ver- 
tex distance. Finally only the decays happening between 
1050 ns and 9000 ns are selected. Earlier tracks might 
overlap with DC signals from the muon. DC signals from 
later tracks may occur after the end of the event record- 
ing. 

It is important to recall that exactly the same algo- 
rithms are applied to data and simulation, reducing the 
dependence of our muon decay results on the precision of 
the algorithms. The evaluation of the systematic uncer- 
tainties from the detector response is accomplished using 
event selection criteria identical to that of the analysis, 
and therefore integrates the effect of the cuts in the un- 
certainties (Sec. |V|) . 



D. Muon decay parameter fit 

The p-9 spectra obtained for data and simulation are 
now compared to perform a momentum calibration and 
to extract their difference in terms of decay parame- 
ters. The data and simulation spectra have very different 
muon decay parameters, compared to the precision of the 
measurement, because of the hidden parameters in sim- 
ulation. This difference is typically a few parts in 10~ 3 
and it biases the edge fit of the momentum calibration 
performed at the kinematic end-point of the two spectra. 
The shape of the spectrum near the end-point is sensi- 
tive to this difference, so it is necessary to include in the 
simulation the derivatives weighted according to the re- 
sults from a prior decay parameter fit. For this reason 
the two fitting procedures are applied iteratively, starting 
with the decay parameter fit. Only one iteration of the 
momentum calibration is needed to reach convergence. 

The muon decay parameter fit procedure exploits the 
linearity in the decay parameters r/, p and the products 
P M £ and P^S [Eq. <gj|]. The difference between the data 
spectrum (Sd) and the Monte Carlo simulation spectrum 
(Smc) can be expressed in terms of derivative spectra of 
the decay parameters 17fl. Schematically: 



S 



MC 



dS_ 
dp~ 

dS 



dPut; 



Ap 



dS 
dP^S 



A(P„#), (8) 



where the Aa, (a — p, P^, Pfi^S) are the free parameters 
of the fit. The effect of the detector response on the p-9 
derivative spectra is simulated using the same code as is 
used for the muon decay spectrum. However, unlike the 
decay spectrum, the derivatives are not positive definite, 
and additional sign information must be passed to the 
fitting software. The radiative corrections are already 
taken into account in the simulation. 



Fiducial regions in the p-9 spectrum are defined to 
reduce bias while maximizing resolution and sensitivi- 
ties to the decay parameters. Only the bins whose cen- 
ter is contained in the fiducial regions are used in the 
decay parameters fit. The maximum momentum cut 
(p max = 52.0MeV/c) avoided the region of the spec- 
trum that was used in a momentum calibration proce- 
dure (described below). The longitudinal momentum cut 
(|p™ ln | = 14.0MeV/c) avoided the region where the he- 
lix wavelength was difficult to determine. The require- 
ment I cos 9\ < 0.96 removed small angle tracks where the 
wavelength was poorly resolved, and | cos#| > 0.54 elim- 
inated large angle tracks with less reliable reconstruction 
due to multiple Coulomb scattering as the path length 
through the chambers became too large. The maximum 
transverse momentum cut (p™ ax = 38.0MeV/c) retained 
only the positrons within the instrumented regions of 
the detector. The minimum transverse momentum cut 
{p? ln = lO.OMeV/c) removed tracks where the helix ra- 
dius became comparable to the wire spacing. The up- 
stream and downstream fiducial regions are symmetric 
about cos 9 = (Fig. We studied the stability of 
the decay parameters with respect to the definition of 
the regions by varying by a few percent all the fiducial 
boundaries. These boundaries were slightly modified for 
this analysis compared to the ones used for the interme- 
diate measurement 81. 
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FIG. 3. (color online) Residuals normalized by the statisti- 
cal uncertainty from the muon decay parameters fit between 
simulation and data. Only bins with their center contained 
in the fiducial regions are used in the fitting procedure. 

A x 2 minimization, using MINUIT ^25] , of Eq. © to 
the data is used to determine the muon decay parameter 
differences. The correlations between the parameters as 
returned by the fitting algorithm are 0.19 for p-8, 0.21 
for p-£ and —0.72 for 8-£. The parameter 77 is not part 
of the fit in this analysis because it is strongly correlated 
to the parameters p and P^ (Sec. |T]) and the fiducial 




.V . 





TABLE III. Mean values of the momentum calibration pa- 
rameters with statistical uncertainties. 



q 0.000 



FIG. 4. Measurement of the angle-dependent momentum mis- 
match at the decay positron kinematic end point for set 84, 
taken with the Al target under nominal conditions. 



regions exclude the low momentum part of the spectrum, 
which is the most sensitive to this parameter. The final 
determination of 6 from A(P^S) is only possible using 
the hidden parameters 8h and £/j in the formula: 



6 h £ h + A(P^S) 



(9) 



However, before unblinding it is sufficient to use the SM 
values to estimate AS. The final value of 5 is recalculated 
after the hidden parameters have been revealed. 



E. Momentum calibration 

The momentum calibration exploits the kinematic end- 
point of the decay positron momentum at 52.83 MeV/c 
to measure the mismatch between the data and simu- 
lation detector responses. Because the planar geome- 
try of the TWIST detector, the momentum loss of the 
positrons exiting the target will have a l/cos0 depen- 
dence. Histograms of the edge region with 10 keV/c mo- 
mentum binning and bins in l/cos0 of width 0.0636 in 
the range 0.5 < |cos0| < 0.9 (1.11 < |l/cos0| < 2.00) 
are produced. For each l/cos0 slice the simulated edge 
histogram is shifted in 10 keV/c steps with respect to the 
data histogram. At each step a x 2 statistic is calculated 
using the difference in bin contents between the spectra. 
The resulting x 2 distribution is fitted with a second-order 
polynomial to determine the momentum shift required 
to minimize the x 2 ■ The momentum mismatch between 
data and simulation versus l/cos0 (see Fig. gj is fitted 
independently upstream and downstream with straight 
lines, 



Ap = at/ |cos(0)| -bi] i= (up,dn). 



(10) 



Target 






O-dn 


bdn 




keV/c 


keV/c 


keV/c 


keV/c 


Ag 


1.8 ± 0.5 


-10.0 ± 0.8 


-3.1 ± 1.3 


-1.7 ± 2.0 


Al 


4.8 ± 0.6 


-6.9 ± 0.9 


-0.2 ± 1.4 


-11. 0± 2.3 



A new data p-9 spectrum is produced by applying the 
momentum calibration for each set on an event-by-event 



basis, and the statistical uncertainties and correlations of 
the calibration parameters are propagated to the muon 
decay parameter error budget. Table UHl shows the mean 
values of the momentum calibration parameters. 

The model used for the propagation of the momentum 
mismatch to the entire spectrum depends on the source 
or sources of the mismatch, which could not be uniquely 
identified. For this reason the final muon decay param- 
eter results are the average of the analyses calibrated 
using a shift that was either constant or scaled with mo- 
mentum. Systematic uncertainties associated with the 
momentum calibration are discussed in Sec. IV CI 



F. Drift chamber calibration 

Improvements to the DC calibration procedures have 
been crucial to reach our final precision for the decay pa- 
rameters. First of all the wire time offsets, which correct 
for the different propagation times of the signals from 
different sense wires, were measured directly from the 
decay positrons in the physics data. Previously the wire 
time offsets were determined from special pion data taken 
only at the beginning and the end of run periods, lead- 
ing to a dominant systematic uncertainty from the time 
dependence of these offsets. For this measurement, a 
downstream scintillator was used in addition to the ex- 
isting upstream scintillator. Both scintillators recorded 
the arrival time of the decay positron as a reference. The 
upstream scintillator is an annular shaped positron scin- 
tillator installed around the main muon trigger scintil- 
lator. The downstream scintillator on the other hand is 
installed outside of the steel yoke and covers most of the 
yoke downstream opening. 

The wire time offsets were extracted from the decay 
positron signals after a time of flight correction. The 
algorithm for fitting these time distributions, which are 
broadened by the drift times of the electrons in the DCs, 
was significantly improved. The mismatch of the offsets 
between data and simulation was estimated to be less 
than 0.5 ns channel-by-channel based on the difference 
in the fit parameter describing the steepness of the DC 
signal rising edge. 

The relative misalignment of the DCs was measured 
and corrected in the analysis to improve the reconstruc- 
tion resolution. A special set of data was taken with 120 
MeV/c pions and with no magnetic field. The straight 
tracks produced by the pions traveling through the entire 
chamber stack were reconstructed. At each wire cham- 
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ber, the residuals were used in an iterative process to de- 
termine the misalignment in translation in the direction 
measured by the chamber (perpendicular to the wires) 
and in rotation around the detector axis with a precision 
respectively of 10 /an and 0.03 mrad. The target PCs 
misalignment was also corrected due to their importance 
to measure the muon-positron vertex distance. The mis- 
alignment between the spectrometer and the magnetic 
field axis was measured on muon decay data using a spe- 
cial helix fitting algorithm allowing for a rotated helix 
axis. This measurement was performed 3 times, each 
time that the spectrometer was removed from inside the 
coils. The three misalignments showed remarkable repro- 
ducibility, being consistent within the 0.03 mrad uncer- 
tainty, with an average value of 0.31 mrad in x and 1.15 
mrad in y. 

The previous TWIST analyses used STRs extracted 
from Garfield simulations. This analysis measured ef- 
fective STRs independently for both simulated and real 
data from the time residuals of the helix fitter on decay 
positron tracks 26] . An iterative procedure modified the 
STRs to reduce the time residuals in subcells of the drift 
cell surrounding the sense wire. All the drift cells are av- 
eraged for each plane. The main advantages of the new 
procedure are to correct for a bias from the helix fitter, 
which systematically defines the closest distance of ap- 
proach of the track to the wire to be less than the actual 
ion cluster distance to the wire, and to allow data and 
simulation to be treated in a more equivalent way in the 
analysis. Furthermore the STRs were measured for each 
plane in data to take into account imperfections in the 
DCs construction such as the cathode foil position rela- 
tive to the wires. On the other hand, one set of STRs, 
measured from the simulation, was applied to all the DCs 
in the simulation analysis since in that case the geometry 
is identical for all chambers. 

The position resolution used during the helix fitting 
to weight the residuals was changed from a constant 100 
/im to an ad hoc expression determined by optimizing 
the momentum bias and resolution in the simulation, 

cr(x) = {120 + 5 [sinh(100a: 2 )] } /rai, (11) 

where x is the distance between the wire and the ioniza- 
tion in cm. Equation (|11[) assigns a larger uncertainty to 
hits that are far from the wire, which are affected more by 
diffusion. For x < 0.1cm there is little sensitivity to the 
position resolution function since a left-right ambiguit}0 
dominates. The improved resolution dependence modi- 
fied the weights used for the track fitting and resulted 
in a difference between data and simulation momentum 
resolutions of < 2 keV/c at the kinematic end-point. 



2 Only the drift time and therefore the distance to the wire are 
known. In those conditions, the left-right ambiguity corresponds 
to the difficulty for the reconstruction algorithm to determine on 
which side of the sense wire the track of the particle occurs. 
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FIG. 5. (color online) Integrated AO distributions for the 
silver target module in data and simulation. 



IV. VALIDATIONS 

Many low-level histograms, such as distributions of 
chamber hits and track lengths, were examined to en- 
sure that the simulation accurately reproduced the data. 
Very little tuning of the simulation was required. As men- 
tioned above, none of the physics processes were tuned 
from their GEANT defaults. Because the systematic un- 
certainty for positron interactions was a leading term for 
our intermediate results, this section includes a detailed 
description of the results of special data taken to test the 
ability of the simulation to reproduce positron interac- 
tions in the detector. These data also allow a precision 
test of reconstruction inefficiencies in data and simula- 
tion. A third subsection describes time spectrum fits, 
which tested the purity of the events in the fiducial re- 
gion as positrons from muon decay. 



A. Positron interactions 

A special data set where the muons are stopped far 
upstream in the muon counter and upstream PCs before 
reaching the DCs is used to validate the relevant positron 
interactions in the simulation, independent of the muon 
decay parameters. In this configuration, a positron from 
a muon decay traverses the entire detector and provides 
two track segments, one on each side of the target. The 
comparison is restricted to single upstream and down- 
stream tracks with hits on at least 16 DC planes and on 
an outer PC plane. The positron track is also required 
to pass within 4 cm of the target center, which limits the 
p-9 phase space over which this comparison can be made. 
The fitted tracks return the position and momentum at 
the drift chamber nearest to the stopping target. The 
difference in angle between the two reconstructed tracks 
provides a test of the ability of the simulation to repro- 
duce multiple scattering through the target module. The 
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distribution of the change in angle is presented for the sil- 
ver target module in Fig. :5j The central width and most 
probable value (MPV) of this distribution are obtained 
from a fit to a Gaussian function. To minimize the ef- 
fects of non-Gaussian tails to this value, the fit region is 
restricted to ±ler about its central value. The agreement 
in both width and MPV is shown in Table IIVI 

The second measurement comes from the change in 
momentum, which tests the validity of the simulation's 
positron momentum loss. The measured momentum dif- 
ference shows a 1 /cos 9 dependence due to the planar ge- 
ometry of the detector. The momentum loss is therefore 
studied using the quantity Ap\ cos9\ as shown for the sil- 
ver target module in Fig. [5] Again a truncated Gaussian 
fit is used to determine the MPV and width of the cen- 
tral peak, which measures soft momentum loss processes. 
For Al there is agreement at the 1 keV/c level. For Ag 
there is a 3 keV/c difference in the MPV momentum loss, 
which is within the uncertainty of 3.5 keV/c for the simu- 
lated positron momentum loss [2?| . The high momentum 
loss tail extending 10 MeV/c above the peak and 3 orders 
of magnitude below the peak height is due primarily to 
bremsstrahlung processes. The integrated counts in this 
tail validate the bremsstrahlung rate at the 1% level, in 
agreement with a separate evaluation based on broken 
tracks (Sec. |VB|) . 
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FIG. 6. (color online) Integrated momentum loss (Ap)|cos#| 
distributions for the silver target module in data and simu- 
lation. The right panel shows the distribution of high mo- 
mentum loss events, due primarily to bremsstrahlung. The 
discontinuity between the two panels is because of the change 
in bin size. 



B. Reconstruction Inefficiencies 

The analysis of the far upstream stops data also deter- 
mines the probability of not finding a track in one-half 
of the detector when it is successfully reconstructed in 
the other half. This inefficiency includes the possibility 
that the tracks physically scatter into or out of the fidu- 



cial region but it is dominated by the possibility of not 
reconstructing an existing track. The double difference 
between upstream and downstream halves of the detector 
and between data and simulation would affect the muon 
decay parameters, P M £ in particular. 

A weighted average of the track inefficiency was com- 
piled from the events that fall within the fiducial region 
for both data and the simulation for each target module 
(Table W} . The weighting was defined using the Bayesian 
interval for the ratio of the failed tracks over the total 
tracks for a given bin. Beam positron tracks are localized 
at the cos 9 = 0.94 fiducial boundary and are therefore 
rather sensitive to inscattering and outscattering; for this 
reason they were removed from the calculation. Table IVl 
shows a clear difference in the upstream and downstream 
inefficiencies due to positron interactions in the target, 
but is reproduced by the simulation at the 0.5 x 10~ 4 
level. Positron interactions in the target module or first 
downstream chambers, including annihilation-in-flight, 
large angle scattering, or production of secondaries that 
confound the reconstruction, will produce such a differ- 
ence in our inefficiency measurements. 



C. Time Spectrum Fits 

To check the consistency of data and simulation, of 
time calibration, and the absence of time-independent 
backgrounds in the data, fits of the selected events to the 
time dependence were performed for a typical data set 
and also for a simulation set. The fits included an over- 
all normalization, the degree of initial muon polarization, 
and also a small time-dependent relaxation of the asym- 
metry [l(| • The fit range was from 2 /is to 9 /us following 
muon arrival to avoid a small decay time distribution 
bias below 1 fis from the algorithm that rejected beam 
positron pileup. Assuming zero uniform background and 
the accepted value of the muon lifetime, acceptable fit 
qualities were obtained for events in the decay param- 
eter fit region. The confidence levels are 75% for set 
84 and 6% for the corresponding simulation, using only 
statistical uncertainties. These results confirm that the 
tracks selected by the fiducial region in data are consis- 
tent with a pure sample of positrons from muon decays. 
However, no systematic evaluation of the lifetime mea- 
surement was attempted, as it was beyond the scope of 
our physics goals and not intrinsically relevant to the 
measurement of decay parameters. 



V. BLIND ANALYSIS UNCERTAINTIES AND 
CORRECTIONS 

Most of the systematic uncertainties originate from a 
mismatch in the apparatus or in physics processes be- 
tween the simulation and the experiment. These uncer- 
tainties are evaluated by purposely exaggerating the mis- 
match in a simulation and measuring the change in de- 
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TABLE IV. Properties of integrated momentum loss (Ap|cos#|) and scattering (A9) distributions. The peak and width of the 
distributions were determined using a truncated Gaussian to remain independent of the long, asymmetric tails. Only statistical 
uncertainties are quoted. 

Silver Aluminium 
Ap|cos6>| A6 Ap|cos6»| A6 

Peak keV/c Width keV/c Peak mrad Width mrad Peak keV/c Width keV/c Peak mrad Width mrad 
Data 40.37 ± 0.46 55.46 ± 0.20 -0.00 ± 0.14 21.09 ± 0.08 32.25 ± 0.42 53.28 ± 0.26 0.13 ± 0.15 11.43 ± 0.06 

Simulation 43.36 ± 0.43 54.84 ± 0.26 -0.20 ± 0.11 20.65 ± 0.10 32.98 ± 0.57 52.21 ± 0.25 -0.09 ± 0.12 11.30 ± 0.05 



TABLE V. Weighted average track inefficiencies upstream 
(US) and downstream (DS) of the target within the fiducial 
region used for the decay parameter fit. Only statistical un- 
certainties are quoted. 

Target Detector Inefficiency (xl0~ 4 ) 

half Simulation Data Difference 



Al 

Ag 



US 3.96 ± 0.16 3.74 ± 0.16 0.36 ± 0.23 

DS 5.71 ± 0.18 6.15 ± 0.19 -0.30 ± 0.28 

US 4.54 ± 0.16 3.74 ± 0.11 -0.30 ± 0.20 

DS 7.13 ± 0.18 7.47 ± 0.15 -0.58 ± 0.25 



cay parameters between this modified simulation and a 
nominal simulation. The difference is the sensitivity of 
the decay parameter to this mismatch. The exaggera- 
tion produces statistically well determined sensitivities. 
A factor corresponding to the ratio between the exagger- 
ated mismatch and the estimation of the real mismatch 
is used to rescale the sensitivity to obtain the system- 
atic uncertainty. The sensitivity to a component of the 
analysis can be obtained by comparing via a decay pa- 
rameter fit the spectra from a standard analysis and from 
an analysis with that component exaggerated, using the 
same data for both analyses. This approach, when pos- 
sible, reduces the statistical uncertainties from the sen- 
sitivity evaluation. It relies on the assumption of linear- 
ity of the systematic uncertainties evaluated, which was 
verified to be valid for large uncertainties such as the 
bremsstrahlung production rate (Sec. IV Bp . Special at- 
tention was also given to avoid the double counting of a 
systematic uncertainty as in the case of the momentum 
resolution during the evaluation of the DC STRs (Sec. 
IV Dp . Table |VT1 summarizes the systematic uncertainties 
by categories that typically contain multiple independent 
uncertainties. 

The weighted statistical uncertainties in Table I VII are 
computed from the statistical errors for the two targets, 
weighted according to the target dependent systematic 
errors. The weighted systematic uncertainties are the 
quadrature sum of the target independent systematic un- 
certainties and the appropriately weighted target depen- 
dent systematic uncertainties. 

As described above, we calculate the systematic uncer- 
tainties for p and S due to a mismatch s as (dp/ds)a s and 
(dS/ds)a s , where a s is our estimate of the possible size 



of the mismatch. a s is common to the p and S system- 
atics, so (dp/ds)(dS/ds)(Ts represents a contribution to 
the correlation between p and 5. The correlation for the 
Ag (Al) target measurement is given by the sum of the 
Ag (Al) and target independent correlations normalized 
by the quadratic sum of the Ag (Al) and target indepen- 
dent systematic uncertainties. The final total correlation 
is the sum of the Ag and Al target correlations weighted 
by the statistical weights used to determine the final de- 
cay parameter measurement. 



TABLE VI. Systematic and statistical uncertainties for the p 
and 5 decay parameters. Most of the categories shown here 
are combinations of several independent uncertainties. 



Category 


Uncertainty ( x 10 4 ) 

P s 


Target independent 






Radiative corrections and rj 


1.3 


0.6 


Momentum calibration 


1.2 


1.2 


Chamber response 


1.0 


1.8 


Resolution 


0.6 


0.7 


Positron interactions^ 


0.5 


0.1 


Others 


0.3 


0.4 


Ag target 






Bremsstrahlung rate 


1.8 


1.6 


Stopping position 


2.0 


6.0 


Target thickness 


3.2 


2.2 


Statistical 


1.2 


2.1 


Al target 






Bremsstrahlung rate 


0.7 


0.7 


Stopping position 


0.2 


0.8 


Statistical 


1.3 


2.4 


Weighted systematic uncertainty 


2.3 


2.7 


Weighted statistical uncertainty 


1.2 


2.1 


Total uncertainty 


2.6 


3.4 



a excluding bremsstrahlung 



A. Target independent uncertainties 

Two systematic uncertainties are external to the 
TWIST measurement. The uncertainty on the radiative 
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corrections is given by the effect of the missing leading 
term 0(a 2 ) on the decay parameters (28[. A numerical 
integration of this term in the TWIST fiducial regions 
showed that it has a similar shape, 5 times smaller than 
the 0(a 2 L) term. The spectrum shape of the 0(a 2 L) 
term is used to evaluate the change in decay parame- 
ters which gave a systematic uncertainty of ±0.16 x 10~ 4 
(±0.63 x 10 -4 ) for p (5). The second external uncertainty 
is due to the significant correlation factor of 0.94 between 
the p and the r\ parameters. The impact of this correla- 
tion on the decay parameters is evaluated by performing 
the decay spectra fit with r\ fixed at the world average 
value lowered or raised by 1 standard deviation. The 
changes in decay parameters are used as systematic un- 
certainties and are equal to ±1.05 x 10~ 4 (±0.12 x 10~ 4 ) 
for p (S). 

The chamber response category contains the system- 
atic uncertainties for the STRs and the cathode foil po- 
sition presented below, and also the asymmetry between 
upstream and downstream efficiency, the crosstalk, and 
the wire time offsets. The upstream-downstream asym- 
metry uncertainty is measured by scaling the upstream 
half of the p-9 spectrum with respect to the downstream 
half according to the difference in inefficiencies between 
data and simulation extracted from the far upstream 
stops data (Sec. IIVB[) . The corresponding systematic 
uncertainty for p (5) is ±0.20 x 10~ 4 (±0.75 x 10~ 4 ). All 
crosstalk in the electronics of nearby wires in the drift 
chambers should be removed by the analysis software. 
An upper limit on a potential systematic uncertainty 
due to remaining crosstalk is obtained by disabling the 
crosstalk removal and using the full change of 0.50 x 10~ 4 
(0.10 x 10 -4 ) for p (6) as the uncertainty. The wire time 
offsets are measured using different scintillators for the 
upstream and downstream halves of the detector which 
can lead to an asymmetry. The potential difference in 
this asymmetry between data and simulation (which was 
also calibrated) is responsible for a systematic uncer- 
tainty of ±0.09 x 10~ 4 (±0.44 x 10~ 4 ) for p (S). 

The spectrometer's reconstruction resolution in angle 
and momentum is obtained from the far upstream stops 
data (Sec. I IV I) . The p-9 spectrum is smeared on an event- 
by-event basis to exaggerate the effect of a resolution 
mismatch. The rescaled sensitivity provides a systematic 
uncertainty for the momentum resolution of ±0.56 x 10 -4 
(±0.70 x 10 -4 ) for p (8). The angle resolution mismatch 
leads to a systematic uncertainty ^ 0.1 x 10~ 4 for both 
parameters. 

The positron interaction category (Table IVI|) includes 
the systematic uncertainty for the backscattering of de- 
cay positrons from outside material that adds confusion 
to the track reconstruction. The rate of backscattering 
positrons normalized to the muons stopping in the tar- 
get is used to measure the mismatch in outside mate- 
rial between data and simulation. The systematic effect 
of the outside material is determined by comparing a 
nominal simulation and the simulation matching set 83 
in which the downstream beam package is added (Table 



HI)) . The corresponding systematic uncertainty for p (5) 
is ±0.48 x 10~ 4 (±0.13 x 10~ 4 ). The decay parameter 
difference between sets 83 and 84, respectively with and 
without downstream beam package, is consistent with 
this uncertainty. 

The category of systematic uncertainties under the 
name "others" in Table I VII contains uncertainties that 
are < 0.3 x 10~ 4 for both p and <5. The overall spac- 
ing in z of the wire chamber planes was established to 
a fractional accuracy of 5 x 10" 5 . An analysis with the 
z positions exaggerated by a fractional change of 10 -3 
showed a corresponding change in the momentum cal- 
ibration. The changes in p and 6 were negligible. A 
small correction to the magnetic field is obtained by fit- 
ting an analytic function to the difference between the 
measured field map and the OPERA-3D map. A com- 
parison of set 84 analyzed with this correction with the 
nominal analysis is used to obtain corrections and as- 
sociated uncertainties. This analysis shows significant 
corrections to the energy calibration parameters. How- 
ever, after the new calibration is applied, the change to 
p and S is < 0.1 x 10~ 4 . Finally, the uncertainties for the 
muon and positron beam intensities are also part of this 
category and are negligible. 



B. Bremsstrahlung and <5-electron production rate 

A difference between data and simulation in the rates 
for the emission of bremsstrahlung photons or S electrons 
would affect the decay parameter measurement. Primar- 
ily these processes modify the positron momentum and 
angle between the muon decay vertex and the beginning 
of the tracking region, thus altering the reconstructed p-9 
spectrum shape. Additionally a large change in positron 
momentum within the tracking region can lead to the 
identification of two separate track segments by the re- 
construction algorithm. This second effect reduces the 
reconstruction resolution by shortening the primary de- 
cay positron track, but it can also be used to compare 
the bremsstrahlung and 5-electron production rates for 
data and simulation. 

The bremsstrahlung production rate is evaluated by 
counting the number of events containing two recon- 
structed tracks from a single decay positron. The data 
and simulation counts are normalized to the number of 
muons stopping in the target. The momentum of the 
bremsstrahlung photon is deduced from the momentum 
difference (Ap) between the two tracks and is shown in 
the left panel of Fig. [7] The agreement between data 
and simulation is excellent except near Ap = 0. The 
discrepancy there could be due to the loss of hits from 
corners of the drift cells, which happens more in data 
than in the simulation. These additional hits lead to a 
higher rate of broken tracks with very little momentum 
difference between the two track segments in the sim- 
ulation. Events with a Ap between 15 and 35 MeV/c 
are used for the comparison. The average ratio of the 
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bremsstrahhmg production rates from all the data sets to 
their corresponding simulations is equal to 1.024 ±0.004. 
Although this ratio is measured for the relatively low Z 
materials of the chambers, it is assumed to be applica- 
ble for the full range of materials in the detector. This 
assumption is supported by the target energy loss mea- 
surements (Sec. I IV I) . The bremsstrahlung rate is strongly 
dependent on the target material. Thus the sensitivi- 
ties to the production rate are measured separately for 
each target, from the difference in decay parameters be- 
tween a nominal simulation and a simulation with the 
bremsstrahlung production rate exaggerated by a factor 
of 3. Sec the right panel of Fig. [7] The systematic uncer- 
tainties (Table IVI|) are given by the sensitivities rescaled 
by the factor (3 - 1)/(1.024 - 1) = 83.3. A simulation 
with a smaller exaggeration factor of 2 was also generated 
and analyzed, and its results confirm the assumption of 
linearity of the systematic uncertainty. 




10 20 30 10 20 30 

Momentum difference (MeV/c) Momentum difference (MeV/c) 



FIG. 7. (color online) Number of events with two tracks ver- 
sus the momentum difference between the two reconstructed 
tracks. The left-hand side shows nominal data set 74 and the 
corresponding simulation. The right side shows nominal and 
exaggerated simulations. The bottom plots correspond to the 
ratio of the two distributions above. 

Evaluation of the 6 electron production rate uncer- 
tainty is done similarly to that of the bremsstrahlung. 
The production rate is measured by requiring a third 
track from a negatively charged particle along with the 
two track segments from the decay positron. The mo- 
mentum of the S electron is measured directly from the 
reconstruction of the negatively charged track and used 
to select events with <5 electrons in the momentum range 
(6 < p < 16) MeV/c. The average ratio from all the 
data sets and their corresponding simulations is equal 
to 1.007 ± 0.009. The sensitivities to the 5 electron 
production rate are also evaluated using a simulation 
with a threefold exaggerated rate. The contribution is 
±0.07 x 10" 4 (±0.06 x 10~ 4 ) for p (S) to the Table EH 
positron interaction uncertainties. 



C. Momentum calibration 



1. End points fits 



The momentum mismatch between data and simula- 
tion at the end point, which is assumed to be linear with 
respect to l/cos# based on geometrical considerations, 
is characterized by the parameters a up , b up , cidn, and bdn 
as shown in Eq. (10). However, if one assumes that the 
uncertainties are purely statistical, the linear fits result 
in a total \ 2 °f 212.9 for 168 degrees of freedom, cor- 
responding to a p value from all the data sets equal to 
0.011. Evidently the behavior of the mismatch is not lin- 
ear, possibly due to higher order effects or perhaps some 
underlying fine structure in the momentum spectrum for 
each angle bin. The manifest nonlinearity stems from the 
upstream end-point portions of the fits, while the down- 
stream portions have larger statistical uncertainties that 
could mask any nonlinear behavior. 

To account for this nonlinearity we add in quadrature 
an uncertainty of 1.6 keV/c to the statistical uncertainty 
of the momentum mismatch at each l/cos# bin, in or- 
der to achieve an upstream reduced x 2 of one. When 
propagated to the uncertainties of the decay parameters, 
this results in systematic uncertainties of ±0.58 x 10 -4 
(±0.54 x 10~ 4 ) in p (S). 

The observed offset at the end-point between data and 
simulation is ~10 keV/c. To understand this difference 
quantitatively, a number of sources of systematic correc- 
tions and uncertainties must be considered. Approxi- 
mately 4 keV/c of this total is due to a slightly incor- 
rect scale used for the magnetic field in the simulation. 
Small corrections to the momentum calibration parame- 
ters were calculated from the best values for the target 
thicknesses, magnetic field map, and match of the muon 
stopping distribution (to be described in Sec. I VI Ajl . 
Systematic uncertainties for these parameters have been 
determined using the same simulation studies that were 
used to determine the muon decay parameter uncertain- 
ties. The most significant items are from the magnetic 
field map, the STRs, the z spacing of the chambers, the 
match of the muon stopping distribution, and the target 
thicknesses. Additionally, there is an uncertainty from 
soft momentum loss in the simulation, consisting of an 
ionization momentum loss uncertainty of 2% and a radia- 
tive momentum loss uncertainty of 3% (27j for the drift 
chamber and target materials used. After corrections, 
the magnitude of the mean slope parameters for each 
target is less than 5 keV/c, and the mean offset magni- 
tude is less than 7 keV/c, both with systematic errors 
of « 5 keV/c. This level of agreement shows acceptable 
consistency of the data with the expected accuracy of the 
simulation. 
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2. Propagation model 

The momentum mismatch between data and simula- 
tion is measured only at the kinematic end-point but is 
corrected over the entire spectrum. The predicted mo- 
mentum dependence of this calibration depends on the 
source of the momentum mismatch between data and 
simulation. For instance, a difference in solenoid mag- 
netic field strength leads to a momentum mismatch that 
depends linearly on the momentum and is referred to as 
a scale. Another example is a mismatch in target thick- 
ness, which translates into an angle-dependent shift of 
the momentum (to first order), with the angle depen- 
dence measured by the slopes a up and ddn- Most of the 
observed offset at the end-point could not be attributed 
to a unique source. Therefore it was assumed that the 
propagation of the momentum mismatch is a mixture of 
shift and scale. 

For this reason the decay parameters were computed 
for the two extreme cases of propagation which corre- 
spond to a pure shift with the form 



P corrected P reconstructed 



cost 



or a pure scale, given by 



P corrected 



P reconstructed 



1 



(12) 



(13) 



| cos 6 



where W ell and m 



were defined in the context of Eq. |4] 
The average values of the p and S parameters using the 
shift and the scale propagations are different respectively 
by 2.04 x 10 -4 and 2.16x 10~ 4 . Their mean is used for the 
decay parameter. Half of the difference between shift and 
scale is used as the uncertainty to cover the two extreme 
possibilities. Therefore the systematic uncertainty from 
the propagation model is ±1.02 x 10~ 4 (±1.08 x 10~ 4 ) 
for p (S). 



D. DC STRs 

The accuracy of the helix reconstruction depends on 
the quality of the STRs. In particular, differences be- 
tween the respective accuracies of data and simulation 
STRs can lead to a bias in the decay parameter mea- 
surement. The STRs were derived in both cases using 
the two-dimensional time residual distributions covering 
the entire drift cell (T res ) from the helix fitter. The sen- 
sitivity to a mismatch in STRs is measured by creating 
simulation STRs containing the difference AT res between 
data and simulation STRs. First 44 AT res are created by 
taking the difference between the data and the simula- 
tion T res for each DC. The 44 AT res are fitted with a fifth 
order polynomial function to guarantee the smoothness 
of the STRs created in the next step. Second 44 STR 
tables are created by adding the 44 polynomial functions 



exaggerated by a factor of ten to 44 duplicates of the 
simulation STRs. A simulation is reanalyzed with these 
new STRs and this set is fitted against the unmodified 
set to measure a change in decay parameters. The corre- 
sponding sensitivity of the decay parameters to the STRs 
changes significantly if the propagation model for the mo- 
mentum calibration is a shift or a scale. For this reason 
the sensitivities from both models are averaged and the 
total sensitivity is -7.5 x 10~ 4 (-14.5 x 10" 4 ) for p (S) 
to STRs exaggerated by a factor of 10. 

The momentum resolution at the kinematic end- 
point is very different between the standard and these 
exaggerated-STR analyses of the simulation. However 
the impact of the resolution on the decay parameters is 
already taken into account in a separate systematic un- 
certainty. The systematic effect from the resolution must 
be subtracted from the STR sensitivities evaluated in this 
section to avoid double counting. 

The sensitivities to the reconstruction resolution are 
evaluated from the differences in decay parameters be- 
tween a nominal spectrum and spectra created with the 
events smeared in momentum by different values. This 
procedure is equivalent to a degradation in resolution. 
The contributions of the resolution to the STR sensitiv- 
ities, which must be subtracted from the total sensitivi- 
ties, are — 4.4x 10~ 4 and —4.5 x 10~ 4 for p and S. Finally 
each sensitivity is scaled down by the exaggeration factor 
of 10 to give the systematic uncertainty for the DC STRs 
of ±0.3 x 10~ 4 (±1.0 x 10~ 4 ) for p (S). 



E. Cathode foil position 

The relative position of each cathode foil with respect 
to adjacent anode wires has two effects on the detector 
response. First of all it modifies the electric field and 
consequently the STRs. This effect is included in the 
plane dependent measurement of the STRs in data and 
therefore does not lead to any additional systematic un- 
certainty. The second effect is the change of the drift 
cell size which can change the number of cells crossed by 
each positron. This has an impact on the track recon- 
struction, in particular on the resolution of the left-right 
ambiguity in the helix fitter. 

In the apparatus there are two different sources of un- 
certainty on the foil position. The first uncertainty is 
due to the DC outer foil bulging toward or away from 
the wires as a result of the differential pressure between 
the DCs and the helium-nitrogen gas mixture surround- 
ing the chambers. The permanent foil bulge was toward 
the wires (60±22) pm (average) during the 2006 run pe- 
riod, and away from the wires (8±22) pm during 2007. 
The second source of uncertainty comes from the con- 
struction of the chambers and was estimated to be ±100 
pm on average. 

The sensitivity of the decay parameters to the cathode 
foil position is evaluated by generating a simulation with 
cathode foils moved toward the wires by 500 pm, without 
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modifying the STRs. The fit of this exaggerated simula- 
tion against the corresponding nominal simulation gives 
a sensitivity of 4.0 x 10~ 4 and 5.9 x 10~ 4 for p and 5. In 
this modified simulation the drift cell size is reduced for 
all the planes but in reality some drift cells are potentially 
larger in data than they are in simulation. Therefore the 
average systematic effect is smaller than the estimated 
sensitivity. For this reason the cathode foil position un- 
certainties from the bulge and the chamber construction 
are not added in quadrature but instead only the largest 
uncertainty of ±100 /im is considered, leading to an ex- 
aggeration factor of 5. The corresponding rescaled sys- 
tematic uncertainty is ±0.80 x 10~ 4 (±1.18 x 10~ 4 ) for 
p(S). 

F. Statistics bias correction 

A sensitivity to the difference in statistics between data 
and simulation was discovered in the % 2 minimization 
technique used by the decay parameter fit and the mo- 
mentum calibration. In the situation where the data and 
simulation spectra have the same number of events, the 
difference of the two asymmetric Poisson distributions of 
two bins leads to a symmetric probability distribution for 
the residuals. However all the simulations contain 2 to 3 
times more events than their corresponding data set to 
reduce the statistical uncertainty for the decay param- 
eters. This creates an asymmetric distribution for the 
residuals and a bias in the % 2 minimization [2j| . 

The biases of the decay parameter fit and the momen- 
tum calibration fit were evaluated by performing the fits 
between a data set and subsets of the simulation with 
matching statistics. For each fit parameter, the differ- 
ence between the average of the subsets and the results 
using the whole simulation corresponds to the bias. Cor- 
rections of —0.20 x 10~ 4 and —0.05 x 10~ 4 for p and 
S were applied to account for the average fitting bias of 
the decay parameters. The momentum calibration fitting 
bias corrections were applied differently, on a set by set 
basis, and on the decay parameter measurements from 
the shift and the scale propagation of the calibration to 
the spectrum. The values of the corrections range be- 
tween -0.92 x 10~ 4 and -1.36 x 10~ 4 (-0.31 x 10~ 4 and 
-0.53 x 10~ 4 ) for p (5). 

VI. POSTBLIND ANALYSIS 

The hidden parameters of the blind analysis were re- 
vealed once the differences of the three decay parameters 
between data and simulation were confirmed and the sys- 
tematic uncertainties were fully evaluated. The results 
for the three parameters were consistent with the SM. 
However, the product P^S/p was 1.00192 ±£ooo66- Al- 
though the sign of the deviations of the individual decay 
parameters from the SM is not constrained in the gener- 
alized matrix element treatment the product P^S/p 



must be < 1. This product can be identified with the 
asymmetry between the extremes of cos 6 = ± 1 and eval- 
uated at x — 1 by using Eqs. ([!]), (O, and ©■ 

A measurement of P^S/ p > 1 could have been due to 
the matrix element treatment or the momentum-angle 
functional form being inadequate to describe the data, 
but it could also have been due to a systematic uncer- 
tainty or correction missing or not evaluated properly in 
the analysis. Furthermore, P^S/p from the blind anal- 
ysis was different for the Ag and Al target data by 3.8 
fj. Both the large value of P^S/p and the mismatch 
between Ag and Al target data triggered an exhaustive 
review of the blind analysis and special scrutiny of var- 
ious systematic effects that could explain these results. 
Among the tests performed, effects such as /i + — > e + X° 
decays (where X° is a long-lived unobserved particle) , an 
incorrect value of the rj parameter, or plausible errors in 
the radiative correction implementation, did not resolve 
the mismatch. However we found two corrections that 
were missed during the blind analysis. 



A. Additional systematic uncertainties and 
corrections 

The effect of the muon radiative decay on the p-0 spec- 
trum is included in the radiative corrections of the decay 
positron spectrum of the simulation (Sec. IIII Ap . How- 
ever we neglected to simulate the photon from radiative 
decays. Although the wire chambers are insensitive to 
photons, the electrons and positrons from pair produc- 
tion or Compton scattering of the photons can affect the 
track reconstruction. These processes occur at different 
rates in Ag and Al and therefore potentially bias the Ag 
and Al measurements differently. The effect on the de- 
cay parameters is measured using two simulations of pure 
muon radiative decay using the Fronsdal and Uberall for- 
mula 30] to calculate the momentum and angle of the 
decay positrons and the photons. One simulation con- 
tains all the standard physics processes while the second 
simulation does not include the pair production and the 
Compton scattering of the photon so that, as in the nom- 
inal simulation, the radiative decay photons are absent. 
The decay parameter difference is renormalized using the 
branching ratio of (1.4 ±0.4)% from Q. The corrections 
to p and 5 for the Ag data are 0.59 x 10~ 4 and 0.76 x 10~ 4 , 
and they are negligible for the Al data. 

The second category of corrections and related refine- 
ment of systematic uncertainties is due to the large sensi- 
tivity to energy loss through bremsstrahlung emission by 
the positron as it travels through the target. We expected 
the momentum calibration to correct for a mean muon 
stopping position difference (MSPD) between data and 
simulation. A match with precision at the level of 1 /im 
is required, but the kinematic edge was not very sensitive 
to the large changes in momentum due to bremsstrahlung 
that can affect the spectrum in the fiducial region. Also, 
we assumed that a up = a^n = corresponded to MSPD 
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FIG. 8. (color online) Distributions of number of muon 
tracks ending in a wire plane, versus plane number (top 
panel), normalized to number ending in PC 6 (plane number 
28), for different data and simulation conditions. The bottom 
panel shows the differences of selected pairs of distributions. 



= 0, which turned out not to be true. An improved tech- 
nique was developed to determine MSPD. For the blind 
analysis the systematic uncertainty for a mismatch in 
bremsstrahlung production rate was evaluated only for 
the Ag target and applied to all the data. In the post- 
blind procedure, it is evaluated for each target, and a 
separate uncertainty is added for a mismatch in target 
thickness. 

A measurement based on the distribution of the last 
wire plane hit by the muons is used to evaluate the 
muon's MSPD between data and simulation. The last 
wire plane distributions are normalized to the number of 
muons stopping in the target defined by the counts in 
PC 6 (wire plane 28) which is located just upstream of 
the target (Fig. [5]). The differences between data and 
simulation last wire plane distributions show agreement 
at the percent level for all planes for data sets 68, 74 and 
76, which confirms a match in the mean and the widths 
of the stopping distributions. However, disagreements 
for other sets identified a sensitivity to MSPD. Different 
methods were tested to establish MSPD with improved 
precision. It was found that MSPD could be measured 
using an average of the PC 5 and PC 7 fractional dif- 
ferences (wire planes 27 and 29) where the sensitivity is 
the highest. We verified that MSPD measurements from 
other planes are consistent with the measurement from 
PC 5 and PC 7. The relationship between the average 
fractional difference and the MSPD is extracted from the 
comparison of simulations with known non-zero MSPDs 
such as between the simulations of the data sets 68 and 
74 shown Fig. [5J Each data set and its corresponding 
simulation were compared and MSPDs of up to 1.6 pm 



were determined for the Ag target and 3.8 pm for the Al 
target. The sensitivity of the decay parameters to MSPD 
is determined by creating p-9 spectra for different depth 
intervals in the target using the true stopping position 
of the muons in the simulation. Set by set corrections 
are applied and range from 0.0 to —3.3 x 10~ 4 (0.0 to 
-9.8 x 10~ 4 ) for p (5). Although MSPD is larger in Al, 
the largest corrections are for the Ag target (set 75) be- 
cause of a higher density and bremsstrahlung production 
rate in Ag. We estimate the MSPD uncertainty to be 
1 pm for Ag and 2 pm for Al. The systematic uncer- 
tainties on the correction, determined from the sensitiv- 
ity to MSPD, are respectively for the Ag and Al targets 
±2.0 x 10~ 4 (±6.0 x 10~ 4 ) and ±0.2 x 10~ 4 (±0.8 x 10~ 4 ) 
for p (5) (under "stopping position" in Table fVTl) . 

The accurate measurement of target thickness is based 
on a destructive test that could only be performed after 
the experiment. It gives a thickness of (30.9±0.6) /im and 
(71.6 ± 0.5) pm for the Ag and Al targets respectively. 
The simulation used the prior estimate of the Ag (Al) 
target thickness of 29.5 (71.0 pm), which was based 
on measurements of material samples that were similar to 
the targets, but obviously not identical. The impact on 
the decay parameters of the mismatch in Ag target thick- 
ness was determined by generating a simulation with a 
65 pm thick target. The corresponding systematic un- 
certainty for p (S) is ±3.2 x 10~ 4 (±2.2 x 10~ 4 ). On 
the other hand, the systematic uncertainties for the mis- 
match in Al target thickness were negligible with values 
< 0.3 x 10~ 4 . 



B. Results 

The final results are extracted from all data sets iden- 
tified as valid for analysis of p and 6. These sets are un- 
changed from the blind analysis. However, the postblind 
result for P^£ includes correlation information from the 
measurement of 8 in five data sets not used for P££, 
which reduces its statistical uncertainty but does not 
change the central value. This analysis has already been 
published separately [Toj . The consistency is shown in 
Fig. [5] for measurements taken under various experimen- 
tal conditions (Tables HI and HI)) , demonstrating that the 
simulation reproduces these conditions accurately. The 
Ag and Al data are fitted separately and then are com- 
bined using the target dependent systematic uncertain- 
ties. 

The additional corrections and systematic uncertain- 
ties determined during the reevaluation of the analysis 
(Sec. I VI Al) change the central values of p and S by 
— 1.4 x 10 -4 and —2.3 x 10~ 4 ; both changes are less than 
the total assessed systematic uncertainties, which them- 
selves changed by less than 0.6 x 10~ 4 . The modified 
values of P*£5/p for Ag and Al data are now consistent 

within ss la, while P^S/p = 1.00179 ±0.00071 has de- 
creased but remains somewhat greater than unity. The 
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FIG. 9. Results of the p and 6 measurement for each data set, 
fitted separately for the Ag sets (left panels) and for the Al 
sets (right panels). Only statistical uncertainties are shown 
here. 



TABLE VII. Results of a new global analysis of muon decay 
data, including the present measurements (parameter defini- 
tions in [ll|]). Best fit values and 90% confidence limits are 
given. P£ — 1 is assumed. Global analysis values of the decay 
parameters are also listed. 



Parameter 


Global analysis results (xlO 3 ) 


Qrr 


< 0.30 (0.16 ±0.11) 


Qlr 


< 0.63 (0.39 ±0.18) 


Qrl 


< 44 (25 ± 13) 


Qll 


> 955 (974 ± 13) 


Brl 


< 11 (6.5 ±3.3) 


Blr 


< 0.52 (0.30 ±0.15) 


a/ A 


0.1 ± 1.4 


P/A 


1.4 ±2.4 


a' /A 


-0.1 ± 1.4 


p/A 


-0.5 ±2.4 


Parameter 


Global analysis results 


P 


0.749 60 ± 0.00019 


S 


0.749 97 ± 0.000 28 


i 


0.998 97 ±0.000 46 


V 


-0.002 7 ± 0.005 



final TWIST results for p and S are 

p = 0.749 77 ± 0.000 12(stat.) ± 0.000 23(syst.); 

5 = 0.750 49 ± 0.000 21(stat.)± 0.000 27(syst.); 

These results represent an improvement of a factor of 
respectively 14 and 11 over the pre- TWIST direct mea- 
surements. They are consistent with the SM predictions 
of p = d = 0.75 and furthermore agree with previous 
measurements (Fig. HO)) . 
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FIG. 10. Summary of the central values and total uncertain- 
ties of p and 8 from this analysis, along with the previous 
published measurements f3a - l39l | . 



VII. THEORETICAL IMPLICATIONS 

A. Global Analysis of Muon Decay 

A new global analysis of all available muon decay data 
has been performed including the final TWIST results 
for the decay parameters and their correlations All 
other input values are the same as in the analysis of [Tll | . 
The global analysis used a Monte Carlo method similar to 
that of [3lJ to map out the joint probability distributions 
for 10 variables (see Table Ivn)) . each of which is a bilinear 
combination of the weak coupling constants g]^. The 
constraint of Qrr + Qrl ± Qlr + Qll = 1 is applied 
(see Eq. ([5])), resulting in 9 independent variables; the 
best fit values and 90% confidence limits are given in 
Table IVHI The decay parameters could then be written 
in terms of these independent variables, and the results 
are included in Table IVIII The present analysis makes 
significant improvements in the limits on Qrr, Qlr, and 
Blr compared to the 2005 analysis, and tightens several 
of the other limits. 

The results from this global analysis can be used in 
Eq. @ to place limits on the magnitudes of the weak 
coupling constants |<?2/J; the exceptions are IqYl] ano - 
\gf L | , which are determined more sensitively from inverse 
muon decay, e~v ll — > These limits are presented 

in Table IVIIII Tighter limits from the present analysis 
of up to a factor of 2 compared to the 2005 analysis are 
found in \g\ R \ and \g\ R \. 

A new indirect limit on the value of P££,5/ p can be ob- 
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TABLE VIII. Limits on the weak coupling constants. (Limits 
on \9ll \ and \gl L \ are from Ref. 0].) 



\g s RR \ < 0.035 


Iflflfll < 0.017 


T 


= 


\g s LR \ < 0.050 


IffZfll < 0.023 




< 0.015 


\g^ L \ < 0.420 


\gR L \< 0.105 


T 


< 0.105 


\gL\ < 0.550 


| 5 I L | > 0.960 




= 



tained from the global analysis. The linear combination 

P-ZS = ^Q RR + 2(Q LR -B LR ), (14) 

combined with the constraints < Qrr, < B^r < 
Qlr, an d P£ < 1, builds in the physical condition 
P££6/p < 1, which is required to avoid a negative muon 
decay probability near the end point [see Eq. (Q}]. We 
find P^S/p = 0.999 47 ± 0.000 28 or P^S/p > 0.999 09 
(90% C.L.). This is a significant improvement over the 
previous limit of P*£6/p > 0.996 82 

The quantity Q R — Qrr + Qlr. represents the total 
probability for a right-handed muon to decay into any 
type of electron, a process forbidden under the SM weak 
interaction. The new limits on Qrr and Qlr. shown in 
Table IVTT1 yield a new 90% confidence limit upper bound 
on the combined probability Q R < 0.000 82, a factor of 
6 improvement over the limit from the pre- TWIST num- 
bers. 



B. Neutrino mixing 

It is now established that flavor mixing occurs in the 
neutrino states [HI, [34| . Thus the neutrino state sum- 
mations in Eq. @ need to extend over the additional 
kinematically allowed states and mixings since the ma- 
trix elements of Eq. ([T]) are evaluated in the flavor ba- 
sis. Doi et al. [35| have calculated the dependence of 
the decay parameters on the neutrino masses and mix- 
ings when only chiral vector weak couplings are allowed. 
These calculations have shown that, when all of the neu- 
trinos involved are much lighter than the muon mass, as 
is the case for v e , f M , and v T , the decay parameters in 
Eqs. ([5]) and ^ are unaffected. However if there are 
additional mixed neutrino states with large right-handed 
Majorana masses, the decay rate is modified. For see- 
saw model extensions of the SM, the effect on the muon 
decay parameters is below the precision of the present 
measurement. 



C. Nonlocal tensor interaction 

The coupling constants g]^ R and g£ L are set to zero in 
the general 4-fermion interaction [Eq. (JTJ] because their 
corresponding matrix elements cancel out. However by 
abandoning locality [1, @, one can redefine the tensor 



interaction as 

T T ®T T = \o^®a 0X -^4- ( 15 ) 
l q 

where is the momentum transfer of some virtual bo- 
son. This form of the tensor interaction permits nonzero 
contributions from g RR and g^ L , in addition to g RL and 
g\ R . If the tensor interaction couples equally to quarks 
and leptons, pion decay data require g^ L , <?^, and gJ jR 
to be very small [H, [36[ . The seesaw mechanism to gener- 
ate neutrino masses would require these same three cou- 
pling constants to be identically zero in muon decay. For 
these reasons, Chizhov @ explored how the muon decay 
spectrum would be changed if the standard model were 
augmented by the addition of a single additional coupling 
constant, gRR- 

Nonzero g RR requires the introduction of a new muon 
decay parameter, k (sa 3^), in the differential muon 
decay spectrum, such that Eqs. ([5|) and become: 

2 

F is (x) = x(l -x) + -p (Ax 2 - 3x - xl) 
y 

+r]Xo(l-x) + KXa + F^ C {x), (16) 

Fas(x) = \^x 2 -xl [1 - x + P(Ax ~ 3 

+ {\jl-xl - 1))] + nx Q (2 - x) 
HFlUx). (17) 

The k term in Eq. (fT6| introduces a negligible distortion 
in the isotropic distribution relative to the precision of 
TWIST. In contrast, the linear k term in Eq. (|17l) repre- 
sents a significant modification to Fas ■ Chizhov 5] finds 
that it increases the integral forward-backward asymme- 
try by a factor 18a;oK. 

The three decay parameters measured by TWIST also 
receive direct contributions from gRR. 

P=\{l-2n 2 ) 1 £=1 + 2k 2 , 
C5= |(l-4^ 2 ), 5=^(1 -6k 2 ). (18) 

Thus, k provides both linear and quadratic modifications 
to the muon decay spectrum. This combination implies 
the linear fitting procedure used in this analysis [Eq. ©] 
cannot be altered to fit k directly. 

A study [37| was performed using the theoretical p-6 
spectrum to determine how a nonzero value of k would 
distort the values we obtain for p, 5, and P M £. We find 

3 

Peff ~ |(! - 0Ax k - 2k 2 ), 
( p ^)eff « 1 + 16.5x k + 2k 2 , 
(P^S)eff « |(1 - 1.2x k - 4k 2 ). (19) 
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The latter two equations imply 



S eff = 



17.7xqk.-6k 2 ). (20) 



When Eqs. (IT9")) and (|20l) are combined with our mea- 
sured values for p, <5, P^, and their correlations to calcu- 
late the probability distribution for k, we find —0.009 < 
k < +0.000 5 (90% C.L.). 



For <5, the interval has been more than 20 years since the 
last measurement [30|. Our final results supersede our 
intermediate values. They are consistent with SM pre- 
dictions, placing more stringent limits on physics beyond 
the SM in the weak interaction. 
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